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ABSTRACT 

^ In this paper I show how the statistics of the gravitational field is changed when the system is 
' characterized by a non-uniform distribution of particles. I show how the distribution functions 
' W{F) and IF(F,dF/df), giving the joint probability that a test particle is subject to a force F 
and an associated rate of change of F given by dF/dt, are modified by inhomogeneity. Then I 
' calculate the first moment of dF/dt to study the effects of inhomogenity on dynamical friction. 

oo ; 

, Finally I test, by N-Body simulations, that the theoretical W(F) and dF/dt describes correctly 
the experimental data and I find that the stochastic force distribution obtained for the evolved 
, system is in good agreement with theory. Moreover, I find that in an inhomogeneous background 



the friction force is actually enhanced relative to the homogeneous case. 
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1. Introduction 

Several authors have stressed the importance of stochastic forces and in particular dy- 
namical friction in determining the observed properties of clusters of galaxies (White 
1976; Kashlinsky 1986,1987) while others studied the role of dynamical friction in the or- 
bital decay of a satellite moving around a galaxy or in the merging scenario (Bontekoe & 
van Albada 1987; Seguin & Dupraz 1996; Dominguez-Tenreiro & Gomez-Flechoso 1998). 
The study of stochastic forces is not only important in the framcvi^ork of galaxy forma- 
tion picture, but it is also important for the study of particular aspects of the evolution 
of a number of astronomical systems, such as galactic nuclei, cD galaxies in rich galaxy 
clusters. 

Send offprint requests to: A. Del Popolo, E-mail:antonino. delpopolo@boun.edu.tr 



2 A. Del Popolo: On the nature of gravitational forces 

The study of the statistics of the fluctuating gravitational force in infinite homo- 
geneous systems was pioneered by Chandrasekhar & Von Neumann in two classical 
papers (Chandrasekhar & Von Neumann 1942, 1943 hereafter CN43) and in several 
other papers by Chandrasekhar (1941, 1943a, 1943b, 1943c, 1943d, 1943e, 1944a and 
1944b). The analysis of the fluctuating gravitational field, developed by the authors, 
was formulated by means of a statistical treatment. Two distributions are fundamental 
for the description of the fluctuating gravitational field: 

1) VF(F) which gives the probability that a test star is subject to a force F in the 
range F, F+ dF. W(F), known as Holtsmark's law (Holtsmark 1919), in the case of a 
homogeneous distribution of the stars, gives information only on the number of stars 
experiencing a given force but it does not describe some fundamental features of the 
fluctuations in the gravitational field such as the speed of the fluctuations and the 
dynamical friction. 

2) VF(F, f ) which gives the joint probability that the star experiences a force F and a 
rate of change f, where f — dF/dt. The features of gravitational field that the VK(F) 
does not describe are described by this second distribution W{F,i). Hence, for the 
definition of the speed of fluctuations and of the dynamical friction one must determine 
the distribution W^(F,f). 

The function PF(F) for a homogeneous system was obtained for the first time by 
Chandrasekhar & von Neumann (1942) under the hypotheses that the stars are 
distributed with uniform density in a spherical system, that there are no correlations 
and that — constant when — s- oo and i? — s- oo, where N is the total number of 
stars and R the radius of the stellar system. They also showed that the force probability 
distribution is given by Holtsmark's lawQ Spherical symmetry is used so that at the 
center of the system Fpot = and Ftot — F stock, while the absence of correlation is 
required in order to be able to decompose the gravitational field into a mean and a 
stochastic component. 

Chandrasekhar's theory (and in particular his classical formula (see Chandrasekhar 
1943b)) is widely employed to quantify dynamical friction in a variety of situations, 
even if the theory developed is based on the hypothesis that the stars are distributed 
uniformly and it is well known that in stellar systems, the stars are not uniformly 
distributed, (Elson et al. 1987; Wybo & Dejonghe 1996; Zwart et al. 1997) in galactic 
systems as well, the galaxies are not uniformly distributed (Peebles 1980; Bahcall 
& Soneira 1983; Sarazin 1988; Liddle, & Lyth 1993; White et al. 1993; Strauss & 
Willick 1995). It is evident that an analysis of dynamical friction taking account of 
the inhomogeneity of astronomical systems can provide a more realistic representation 



^ Historically Holtsmark's law was calculated to obtain the probability of a given electric field 
strength at a point in a gas composed of ions. 



A. Del Popolo: On the nature of gravitational forces 



3 



of the evolution of these systems. Moreover from a pure theoretical standpoint we 
expect that inhomogeneity affects all the aspects of the fluctuating gravitational field 
(Antonuccio & Colafrancesco 1994). Firstly, the Holtsmark distribution is no longer 
correct for inhomogeneous systems. For these systems, as shown by Kandrup (1980a, 
1980b, 1983), the Holtsmark distribution must be substituted with a generalized form 
of the Holtsmark distribution characterized by a shift of l/l^(F) towards larger forces 
when inhomogeneity increases. This result was already suggested by the numerical 
simulations of Ahmad & Cohen (1973, 1974). Hence when the inhomogeneity increases 
the probability that a test particle experiences a large force increases, secondly, W^(F, f ) 
is changed by inhomogeneity. Consequently, the values of the mean life of a state, the 
first moment of f and the dynamical friction force are changed by inhomogeneity with 
respect to those of homogeneous systems. 

In this paper I show how the M^(F) and l^(F,f) distributions are changed by inho- 
mogeneity and the effect of this change on dynamical friction. 

The plan of the paper is the following: in Sec. 2 I show how inhomogeneity changes 
M^(F) and what is the effect of having a finite number of bodies in the system. In Sec. 3, 
I calculate M^(F, f ) and its first moment in inhomogeneous systems. In Sec. 4, I show the 
effect of inhomogeneity on dynamical friction. In Sec. 5, I check the theoretical results 
with N-body simulations and Sec. 6 is devoted to conclusions. 

2. The distribution function W(F) in inhomogeneous systems 

The force per unit mass acting on a test star of a finite N-body gravitating system is 
given by the usual formula: 



where the sum is extended to the N stars. Mi is the mass of the i-th field star and 
is its distance relative to the origin. The value of F at a given time depends on the 
instantaneous positions of all the other stars and therefore is subject to fiuctuations as 
these position change. Even if the stellar distribution were constant and homogeneous the 
total force experienced by a star, 'Ptot, would fluctuate around an average value due to 
local Poisson fluctuations in the number density of neighboring stars. Under conditions 
typically found in astrophysical systems like globular clusters and clusters of galaxies, 
the average force 'Ptot, acting on a test star can be naturally decomposed as: 



The first term on the right hand side of Eq. [2] is the mean field force produced by 
the smoothed distribution of stars in the system and can be obtained from a potential 




(1) 



^ pot ^~ stoch 



(2) 
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function, <i>(r), which is also connected, by the Poisson equation, to the system mass 
density, p(r). In the case of a truncated system having a power-law density profile: 

pM = Po(^)^ 0<r<R (3) 

the mean internal gravitational force is given by: 

F.o. = -^r^-^ (4) 

where G is the Gravitational constant, Mtot is the total mass of the system and R is its 
radius. The second term on the right hand side of equation Eq. [2] is a stochastic variable 
describing the effects of the fluctuating part of the gravitational field. This component 
arises because of the discreteness of the mass distribution: it would be null in a continuous 
system. The fluctuating part of Ftot that I have previously indicated F stock can be studied 
by probabilistic methods defining a probability distribution of stochastic force W {F stock)- 
In the case of an infinite homogeneous system, under the other hypotheses given in Sect. 
1, the random force distribution law is given by Holtsmark's distribution: 

N3/2- 



W{Fstock) = r tdtsm{tFstoch)exp \-n{Gmt)^/'^^ 

Jo I 15 



(5) 



(Chandrasekhar & von Neumann 1942), where m is the mass of a field star, n the mean 
density and G the gravitational constant. This equation shows that the stochastic force 
probability distribution depends only on the mean density n or, equivalently, on the static 
configuration of the system. 

The stochastic force probability distribution in inhomogeneous systems can be calcu- 
lated in the same way and under the same hypotheses as that in homogeneous systems. 
Kandrup (1980) obtained an equation for the distribution of random force of an infinite 
inhomogeneous system with probability density T(r) = The resulting equation is a 
generalization of the Holtsmark law: 

"/^.„^N(3-p)/2 f°° dz{z-smz) 



WiFstock) = ^^^ r tdtsm{tFstoch)e^p -^(Gmi)(3-p)/2 f 

Jo . ^ Jo 



(6) 



where a = ^^~^l^p)^''' = 47ra (because the total number of stars at a distance R is given 
by: N{R) — in drr^T{r) = 47rai?('^^P'/(3 — p)). Equation [6] reduces to the Holtsmark 
distribution for a uniform system in the case p — 0. The probability density must be 
chosen in such a way to ensure a convergent integral in Eq. [6l this restricts the choice 
to p < 3. If the system is finite, as e.g is the case for globular clusters, the Holtsmark 
and Kandrup laws must be substituted by a law in which N is finite. This last can be 
calculated as follows: suppose to have a cluster of radius R containing N stars of mass m 
with probability distribution law given by T(r) — define the characteristic function 
of the probability distribution M^(F) B 

C(t) = / exp(itF)VF(F)dF (7) 



Note that from now on, F stands for Fstoch 
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Suppose now that the system contains only one star. If W{Fi) is the probabihty distri- 
bution of the force due to the star calculated at the origin one verifies that: 

W{F,)dF, = r(r)dr (8) 

Using the usual equation = ^^r^ I obtain the probability distribution for one star: 



a{GmtY^-P^/^ GM 



< Fi < oo 



(9) 



2|r.|(9-p)/2 ^2 

Introducing Eq. [9] into Eq. [7] the characteristic function for the force due to one star 
turns out to be given by: 



M(t) = Ci(t) 



dz 



where z = |t • F|. Finally the probability distribution for F is: 
1 

and the probability distribution of the modulus of the force is: 



Wn{F) = — 



sm(t F ) . 9 

tF ^ ' 



Wn{F) = — / tdtsinitF) 



dz 



N 



(10) 



(11) 



(12) 



2i?(3-p) 

that reduces to Holtsmark's law for iV ^ oo. In terms of = Gma^/<3-?') ^^^^ 2/ = 
tGma^^ , Eq. [T^can be written as: 



2/3 



ydy sin{Py) 



(3-p)/2 foo 



-, N 



2N 



sm z 

z(7-p)/2 



dz 



(13) 



In the case the system is also clustered the Wn (F) distribution can be calculated following 
Antonuccio-Delogu & Atrio-Barandela (1992) (AA92). I suppose that particles have the 
density distribution: 



•r('') = — exp 



(14) 



where a and p and vq are three constants. The stochastic force distribution is given by: 

dH 

(15) 



Wn{F)= J ^Mjv(t)exp(-*tF) 
The term Mn is given by: 



Miv(t) = 



Anjt) 



^2^ N'A2{t) 



(AA92), where A„(t) is given by: 



An{t) 



N 



(16) 



(17) 



N 2 

A2{t) is given by Eq. [T7]with n = 2 while S(t) is given in the quoted paper (Eq. 34) and 
Integrating Eq. [Tn]I finally obtain: 



N 



27rr^-''r[(3-p)/2] ' 



2F 



Wn{F) = 47r^|F|^I^Ar(F) = — / tsin(tF) 



TT 



1 + 1(1-1)^ 

2^ N>A2{t) 



-(GmO^'-P^/' 



Gmt 



Gmt sin z 
exp( —) 



-I N 



rlz ^z(7-p)/2 



(18) 
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3. The distribution function W(F,f) in inhomogeneous systems 

To calculate W{F,f) in an inhomogeneous system I consider a particle moving with a 
velocity v, subject to a force, per unit mass, given by Eq. ([T]) and to a rate of change 
given by 

AT 

V, 3r,(r,V,) 



dF 



N 



(19) 



where is the velocity of the field particle relative to the test one. 
The expression of VF(F, f ) is given following Markoff''s method by (CN43): 

W{F, f) = ^ A{k, S) {exp [-*(k* + S*)]} dkd-E 

with A(k, S) given by 

A(k,S) = e^"'^'^'''^) 

being 

C(k, S) = / / / r [1 - exp i (k* + S*)] drdVdM 

Jo J —oo J —oo 

where n is the average number of stars per unit volume while ^> and are given by the 
following relations: 

(23) 



(20) 



(21) 



(22) 



|r|3 



T dF ,^ 

dt 



V 3r(rV) 



(24) 



and r(V, r, M)dV drdM is the probability that a star has velocity in the range V, V+dV, 
r, r + dr and mass in M, M + dM . 
Now I suppose that t is given by: 



V'(j'(M)iun 



(25) 



where a is a constant that can be obtained from the normalization condition for r, j a 
parameter (of dimensions of velocity^^), an arbitrary function, u the velocity of a field 
star. In other words I assume, according to CN43 and Chandrasekhar & von Newmann 
(1942), that the distribution of velocities is spherical, i.e. the distribution function is 
V'(u) = ■;/;(j^(M)|up), but differently from the quoted papers I suppose that the positions 
are not equally likely for stars, that is the stars are inhomogeneously distributed in space. 
Following Chandrasekhar, a lengthy calculation leads to find the function A(k, S): 

A(k,S) 



„[- f (GAffc) V^-S(p) + ^fc^(GM) 



|v|-[Sisin7 - 2S3COS7] -^(p) - ^(GM)^fc^^ 

8 



— sm 



-f(i-')l} 



ur 



i (a + 6 + c)(S? + E^) + ^(2c + d) 



|v|^ • [(asin^ 7 + ccos^ 7)E^ + (bsin'^ 7 + ccos^ 7)E2 
(csin^7 + (icos^7)i;3 + /SiEs sin7C0S7]} ] 



(26) 
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where 



u ■ 



(27) 



u is the distribution of the velocities of the field stars, and v is the velocity of the test 
star, 

sin X 3 sin x 3 cos x 



B{p) = 



r(4-p)/2 .7.(8-p)/2 .7.(6-p)/2 



• dx 



z — sm z 
z(7-p)/2 



24 



• dz 



216 



6 = 



p+1 (p+l)(p-3) (p+l)(p-3)(p-7) 
-72 



(p+l)b-3)(p-7) 
-36 



(p-3)(p-7) 
-2 12 18 



p+1 3 — p 7 — p 
-4 48 



288 



p+1 (p+l)(p-3) (p+l)(p-3)(p-7) 



/ 



24 



12 



144 



(28) 
(29) 
(30) 
(31) 
(32) 
(33) 
(34) 
(35) 



p+1 (p+l)(p-3) 3-p (p-3)(p-7) 

7 is the angle between k and v, S = (Ei,I]2,S3). Eq. introduced into Eq. ([^0]) 
solves the problem of finding the distribution l/l^(F,f) and makes it possible to find the 
moments of f that give information regarding the dynamical friction. Eq. (j26p can be 
written in a more compact form, namely: 

^(k,S) = e-"-^^{\- zffp(k,S) 

+ 6/c^^]-[Q(S)+/ci?(S)]} 



if we define: 



a=|(GM)^ •B(p) 
g = |(GAf)^.|v|.A(p) 



p(k, S) = fc 2" . sin 7 — 2E3 cos 7] 



(36) 

(37) 
(38) 
(39) 



H[-Ki+0]-'4-f(f+')l}-H 



Q(S)= i (a + 6 + c)(E2 + E2) + ^(2c + d) 



(40) 
(41) 
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— (asin^ 7 + ccos^ 7)11^ + (6 sin^ 7 + ccos^ 7)23 
+ (csin^ 7 + dcos^ 7)113 + /E1S3 sin7COS7 



As I stressed in the introduction, the study of the dynamical friction is possible when 
we know the first moment of f . This calculation can be done using the components of f 
ifijjjk)- We have that: 
C^WiF,f)f,d{ 



(42) 



W(F) 

The distribution function VF(F), giving the number of stars subject to a force F, can be 
calculated as follows: 



1 pQC pOC pQO 

W{F) = -— / / / {el-^C^P+^f)]} • A(k,S) dkdSdi 
647r Jq Jq Jq 

Integrating, I find: 



W{F) 



1 



2tt^F 



{e 



[_^2^(3-p)/2] ^ 



} ■ksm{kF)dk 



(43) 



(44) 



This equation gives the generalized Holtsmark distribution obtained by Kandrup (1980a 
- Eq. 4.17) and provides the probability that a star is subject to a force F in a inhomo- 
geneous system. 

In order to calculate the first moment of f wc need only an approximated form for A(k, S) 
(see Chandrasekhar 1943): 



A(k,S) = {e 



2l.(3-p)/2, 



} • [1 - igpik, S)] 



(45) 



Using this last expression for A{k, S) and Eq. dSH]), Eq. ([391), Eq. (gS]), Eq. (jiij) . Eq. (|45l) 
and performing a calculation similar to that by CN43 the first moment of f is given by: 



f = - 
where 

m - 



■A{p)-B{p)^-.- 

7tH{I3)i3 2 



3F-V ^' 
V --T^ • F 



(46) 







6/ 


e 


JO 










JQ 




2 /"^ 



(3-p) 

^(x//3)^- 



(3-p) 

,(a;//3)^- 



cosa; 



sma; 

2;(2-p)/2 ~ 

sinx 



dx 



a;(P-2)/2 



dx 



tt(3 , 



xsin[x)dx 



(47) 
(48) 



The results obtained here for an inhomogeneous system are different [see Eq. (j46|) ]. as 
expected, from those obtained by CN43 for a homogeneous system (CN43 - Eq. 105 or 
Eq. (j55p V At the same time it is very interesting to note that for p — Q (homogeneous 
system) the result here obtained coincides, as is obvious, with the results obtained by 
CN43. In fact, for p = Eq. ^ reduces to: 



f = Gm 

6 









p=0 



3F-V ^ 

V - , ,^ • F 



(49) 
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Defining 

H{/3)^^[ ( e [^T^ 2: sin xda; 



(50) 



we have that 

[m]p=o 



6 sinx 



In this way we can write Eq. as: 



f — Gmn 

3 



3-/„^g(/3)d/3 
P ■ H{(3) 



6 cos X + 2x sin x 



3F-V ^ 

V - , ,^ • F 



dx = 2,-K H{l3)dl3--Kl3H{l3){bl) 




(52) 



this last equation coincides with Eq. (105) of CN43. 



In inhomogeneous systems, Chandrasekhar equation 

can be written, using Eq. (|46|) . as: 

,d|FL . a^GML{B) vF 



(53) 



(54) 



In order to check the validity of the quoted relation (Eq. (|54p ). I have performed numerical 
experiments. This was done by evolving 100000 points (stars) acting under their mutual 



gravitational attraction. From the evolved positions and velocities of the stars, 



/rfF\ 

\ dt I 



computed as a function of velocity and force, similarly to Ahmad & Cohen (1974), and 
then compared with Eq. ((54|) as I shall describe in the following. 

As previously quoted, the results obtained in the present paper for an inhomogeneus 
system are different [see Eq. (|46| ]. as expected, from that obtained by CN43 for a ho- 
mogeneous system (CN43 - Eq. 105). In a inhomogeneous system, in a similar way to 
what happens in a homogeneous system, f depends on v, F and 9 (the angle between v 
and F) while differently from homogeneous systems, f is a function of the inhomogeneity 
parameter p. The dependence of f on p is not only due to the fimctions A{p)^ B{p) and to 
the density parameter a but also to the parameter (3 — |F|/Qh- In fact in inhomogeneous 
systems the normal field Qh is given by Qff GM(aB(p)/2)2/(3-p), clearly dependent 
on p. 



4. Dynamical friction in inhomogeneous systems 

The introduction of the notion of dynamical friction is due to CN43. In the stochastic 
formalism developed by CN43 the dynamical friction is discussed in terms of f : 

3F-V 



f = -^GmnB{(3) 



F 



where 



B{f3) 



3 ■ WjW 



(55) 



(56) 
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As shown by CN43, the origin of dynamical friction is due to the asymmetry in the 
distribution of relative velocities. If a test star moves with velocity v in a spherical 
distribution of field stars, namely (t>(u) then we have that: 



V = u - V = -V (57) 



The asymmetry in the distribution of relative velocities is conserved in the final Eq. (|55p . 
In fact from Eq. ([55)1 we have: 

f = |«.).f (58) 

(CN43). This means that when v F > then ^ > 0; while when v • F < then 
< 0. As a consequence, when F has a positive component in the direction of v, 
|F| increases on average; while if F has a negative component in the direction of v, |F| 
decreases on average. Moreover, the star suffers a greater amount of acceleration in the 
direction — v when v • F < than in the direction +v when v • F > 0. 
In other words the test star suffers, statistically, an equal number of accelerating and 
decelerating impulses. Being the modulus of deceleration larger than that of acceleration 
the star slows down. 

At this point we may show how dynamical friction changes due to inhomogeneity. From 
Eq. (|46p we see that ^ differs fr om that obtained in homogeneous system only for the 
presence of a dependence on the inhomogeneity parameter p. If we divide Eq. (|46p by 
Eq. ([55)1 we obtain: 

V / Horn. 

- ■A{p)-B{p)— ^-^^ 59 

If we consider a homogeneous system, p — 0, the previous equation reduces to: 



= 1 (60) 

Horn. 

In the case of an inhomogeneous system, p 7^ 0, we see that: 

hilh. ^ nP/(3~f)i^(/3(n,p)) (61) 



{ft) 
In the 



where 



6~p 



. -i . (i) - . M,, - i,,,, A . JL^r (8^) 

As I show in Fig. 1, this last equation is an increasing function of p. This means that 
for increasing values of p the star suffers an even greater amount of acceleration in the 
direction — v when v • F < than in the direction +v when v • F > 0, with respect to 
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0.5 



1.5 



Fig. 1. The function F{/3) for several values of the inhomogeneity parameter p; solid line 
p = 0, dashed line p = 0.1; dotted line p = 0.5 

the homogeneous case. This is due to the fact that the difference between the amplitude 
of the decelerating impulses and the accelerating ones is, as in homogeneous systems, 
statistically negative, but now larger, being the scale factor greater. This finally means 
that, for a given value of n, the dynamical friction increases with increasing inhomogeneity 
in the space distribution of stars (it is interesting to note that this effect is fundamentally 
due to the inhomogeneity of the distribution of the stars and not to the density n). In 
other words two systems having the same n will have their stars slowed down differently 
according to the value of p. This is strictly connected to the asymmetric origin of the 
dynamical friction. 

In addition, by increasing n the dynamical friction increases, just like in the homogeneous 
systems, but the increase is larger than the linear increase observed in homogeneous 
systems. 



5. N-body experiments 

To calculate the stochastic force in an inhomogeneous system, I used an initial config- 
uration in which particles were distributed according to a truncated power-law density 
profile: 



0<r<R 



(63) 



12 
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— 1 — J « — I 1 — i — i — L — I — 1 1 I t 1 I I 1 I I L 

I 2 3 4 



5 



Fig. 2. The theoretical probability distribution of stochastic force iir an inhomogeneous 
system for different values of N. The solid line is Kandrup's distribution for an infinite 
inhomogeneous system, the long dashed line is Kandrup's distribution for a finite inho- 
mogeneous system virith N = 50, the dotted line is the same distribution for N = f 00. 
The Kandrup's distribution for a finite inhomogeneous system with N — f 000 is indistin- 
guishable from that for the infinite system. The force is measured in units of Gma^/^^^^-'. 

(see Kandrup 1980). If the velocity distribution is everywhere isotropic then the equation 
relating the configuration space density p(r) to the phase space density f(E) is: 



where U (r) is the potential (normalized to zero at infinity) . Eq. (|64|) may be converted 
into an Abel integral equation and inverted, giving the phase space density: 



(Eddington 1916; Binney & Tremaine 1987). The initial conditions were generated from 
the distribution function that can be obtained from Eq. (j65p assuming a cut-off radius 
i?, = 1, the mass of the system M = 1, ro ~ O.IS^ and G = 1. All the particles had equal 
mass. To have a system whose total mass is contained in a unitary sphere, Eq. (|64|) was 
renormalized and consequently also the potential of the system which is obtained from 
Eq. ([M]) through Poisson's equation. The system of 100000 particles was evolved over 
150 dynamical times using a tree N-body code (Hernquist 1987). During the evolution 
of the system the essential quantities such as position, force, etc., for each test point of 
^ This is the value I used, remember however that the distribution is scale-free 




(64) 




(65) 
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1 1 T 1 1 1 1 1 1 1 T 1 1 1 T \ 1 rZl 




0.2 0.4 0.6 0.8 1 

F (units of GM/R®) 

Fig. 3. Experimental distribution of the stochastic force in an inhomogeneous system. 
The sohd Hne is Kandrup's distribution for a finite inhomogeneous system { p = 0.01) 
of 100000 particles, the histogram is the experimental distribution of stochastic force 
obtained from a system of 100000 particles as described in the text. The force is measured 
in units of ^j^. 

the system was sampled every ^ of a dynamical time (see Ahamad & Cohen 1973 for 
details). The stochastic force, F stock, was calculated observing that at the centre of a 
spherical system we have: 

Ftot — F stock + F„ied — F stock (66) 

because the mean field force, F„ied, is equal to zero. The force was calculated on a 
point at the centre of the system because theoretically Kandrup's distribution gives the 
probability distribution of the stochastic force for a particle at the centre only. When 
points displaced away from the centre are used the stochastic force distribution must be 
calculated as follows: 

a) the stochastic force should be calculated subtracting the mean field force from the 
total force: 

Fstock = Ftot — Fmed (67) 

b) the theoretical distribution must be numerically simulated as done by Ahmad & Cohen 
(1973). 

I F I ■ 

The average of the is a function of velocity, v, force, -F, and the angle between 
them. The test of Eq. (|54p was performed in a similar way to that of Ahmad & Cohen 
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Fig. 4. Experimental distribution of stochastic force in a clustered system. The solid 
Hne is the AA(92) distribution for a clustered system of 100000 particles. The histogram 
is the experimental distribution of stochastic force obtained from a clustered system 
of 100000 particles as described in the text. The dot dashed line and histogram are a 
Holtsmark's distribution and an experimental distribution for an homogeneous system 
of 100000 particles. 

(1974), namely by integrating out two of the variables and examining {^^^) against the 
remaining one (see Ahmad & Cohen 1974). As in Ahmad & Cohen (1974), (^^|^)|f| 
indicates {^^^) after integrating out the angle and velocity, while (^^|^)|v| is (^^|^) after 
integrating out the angle and force. In integrating out the angle, one cannot average over 
the entire range to tt, since that would give zero. Instead, the cosine was averaged from 
to 7r/2, and in the numerical experiments only those particles having a cosine in the 
quoted range were used. In order not to waste the statistics for half the particles the 
same trick of Ahmad & Cohen (1974) was used, namely when the cosine is in the range 
of 7r/2 to TT, the sign of is changed and it is coimted in the same statistics. This 
corresponds to assume that the cosine between v and F is uniformly distributed, which 
is what is found in numerical experiments. The average value of the cosine between 
and 7r/2 is 1/2. 

For a general distribution (|v|) can be calculated as usual: 




(68) 



and {^^^^) can be written, in units of 




GM{\v\), as: 




(69) 
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Fig. 5. The average value of the time rate of change of the magnitude of the force as the 
function of the velocity. The solid line refers to the homogeneous case (Chandrasekhar & 
von Neumann 1943). The dotted and dashed line refers to the cases p = 2.5 and p = 4, 
respectively (see Eq. 26). Crosses represent the experimental points. 



In the particular case of a Maxwellian distribution for velocities: 

i^ = ^^^pi-jM^) (70) 

where p = 2^iy ' that: 

(|v|> = ^ (71) 
we have that: 

that expressed in units of ^-^j^ra^-p GM , then Eq. ([72]) becomes: 

Similarly to Ahamd & Cohen (1974), since to integrate out the force from Eq. ((54l) 
one has a divergent result, I consider only particles up to a certain maximum value of 
the force, P^ax- for example in the case p = 0, f3^aK ~ 10.7, that involves 97% of the 
particles. As observed by Ahmad & Cohen (1974), any cutoff of the force can be used 
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Fig. 6. The average value of the time rate of change of the magnitude of the force as the 
function of the force. The solid line refers to the homogeneous case (Chandrasekhar & 
von Neumann 1943). The dotted and dashed line refers to the cases p = 2.5 and p = 4, 
respectively (see Eq. 19). Crosses represent the experimental points. 

as long as it is taken into account in both the experiment and the analytic evaluation of 
(^)|v|- Defining: 

and 



Jl^'- H{l3)d(3 



(Bi m = (75) 



we find, in units of (i) {Bl{(3))a^ GM , that: 

{^)l.\=A{p)-B{p)^p -v (76) 

The results of calculation and numerical experiments are plotted in Fig. 2-6. 

As Ahmad & Cohen(1973) showed, the stochastic force probability distribution for an 
infinite homogeneous system (Eq. [5|) and that for a finite one (Eq. [12] with p = 0) almost 
coincide for N ~ 1000. In Fig. 2, I show that the same result holds in the inhomogeneous 
case {p = 0.5). The different curves are obtained from Eq. [121 describing the stochastic 
force distribution in a finite inhomogeneous system (long-dashed line) , for increasing value 
of N {N = 50, 100, 1000) and from Eq. [6l which gives the stochastic force distribution in 
an infinite inhomogeneous system (solid line). When N ~ 1000 the two distributions are 
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indistinguishable, meaning that for N ~ 1000 the stochastic force in an inhoniogeneous 
system are equivalently described by the Kandrup's law for an infinite (Eq. [5]) or finite 
system (Eq. I12p. In Fig. 3 I compare Kandrup's distribution for a finite inhomogeneous 
system {p — 0.01) with the histogram of forces obtained from a system of 100000 particles. 
The plot shows a good agreement between the theoretical and numerical distribution. 
In Fig. 4, I study the distribution of stochastic force in a clustered system. The solid 
line is the AA(92) distribution for a clustered system of 100000 particles. The histogram 
is the experimental distribution of stochastic force obtained from a clustered system 
of 100000 particles. The dot dashed line and histogram are a Holtsmark's distribution 
and an experimental distribution for an homogeneous system of 100000 particles. The 
plot shows as expected (Prigogine & Severne 1966; Gilbert 1970; AA92) that there is 
an increase in higher force probability with respect to homogeneous and inhomogeneous 
systems. In Fig. 5, I plot the average value of the time rate of change of the magnitude 
of the force as a function of the velocity. The solid line refers to the homogeneous case 
while the dotted and dashed lines refer to the cases p — 2.5 and p = 4, respectively. 
Crosses represent the experimental points. As shown, experimental points follow a linear 
relationship and there is a good agreement with the theoretical prediction, (Eq. I69p. In 
Fig. 6, I plot the average value of the time rate of change of the magnitude of the force as 
the function of the force. As in the previous figure, the solid line refers to the homogeneous 
case while the dotted and dashed line refers to the cases p = 2.5 and p = 4, respectively. 
In this case, the dependence is no longer linear: it behaves like B{f3) in the homogeneous 
case. A comparison with numerical experiments shows that there is a good agreement 
with the theoretical prediction, (Eq. [76|) . The situation described in Ahmad & Cohen 
(1974), that the experimental data were somewhere in between the theoretical curves for 
the one-particle and the infinite-particle case, is no longer present and the agreement 
is better, now. This is due to the larger number of particles used in the simulations. 
The plots show that Chandrasekhar & von Neumann's theory of dynamical friction in 
gravitational systems gives a good description of experimental data (solid line and data), 
and so does the generalization of the quoted theory to inhomogeneous systems (dotted 
line, dashed line, and respective data). In inhomogeneous systems, Chandrasekhar's result 
which relates the frictional force only to the local properties of the background at the 
position of the object, is no longer true, and friction depends on the global structure of the 
system. This point is in agreement with Maoz (1993), who showed that in inhomogeneous 
media the friction, unlike Chandrasekhar's formula, depends on the global structure of 
the entire mass density field. 
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6. Conclusions 



In this paper I showed how the distributions W(F), W(F,f) and the first moment 
f — dF/dt are changed in an inhomogeneous system. I obtained an expression for W(F) 
in finite inhomogeneous systems, an expression for W(F,f) and one relating f = dF/dt 
to the degree of inhomogeneity in a gravitational system. I showed the implications of 
this result on the dynamical friction in a inhomogeneous system and in particular how 
inhomogeneity acts as an amplifier of the asymmetry effect giving rise to dynamical fric- 
tion. Moreover, I tested by numerical simulations the previous results. For what concerns 
W{F), I showed that Kandrup's (1980) theory, describing the stochastic force probability 
distribution in infinite inhomogeneous systems, and AA92 theory, giving the stochastic 
force probability distribution in weakly clustered systems, describe correctly the observed 
behavior. The results shows that in inhomogeneous and clustered systems there is an in- 
crease in higher force probability with respect to homogeneous systems. Furthermore I 
showed that for N > 1000 Kandrup's theory can be applied to finite systems of particles. 

In agreement with Ahmad & Cohen (1974), the stochastic theory of dynamical friction 
developed by Chandrasekhar & von Neumann (1943), in the case of homogeneous grav- 
itational systems, gives a good description of the results of numerical experiments. The 
stochastic force distribution obtained for inhomogeneous systems, in the present paper, is 
also in good agreement with the results of numerical experiments. Finally, in an inhomo- 
geneous background the friction force is actually enhanced relative to the homogeneous 
case. 

References 

Ahmad A., Cohen L., 1973, ApJ, 179, 885 
Ahmad A., Cohen L., 1974, ApJ, 188, 469 

Antonuccio-Dclogu, V., Atrio-Barandela, P., 1992, Ap.J. 392, 403 
Antonuccio-Delogu V., Colafrancesco S., 1994, ApJ, 427, 72 
Bahcall N.A., Soneira R.M., 1982, ApJ, 262, 419 
Bontekoe, T. R., van Albada, T. S., 1987, MNRAS, 224, 349 

Binney J., Tremaine S., 1987, Galactic Dynamics, in Princeton Series in Astrophysics, Princeton 

University Press. 
Chandrasekar S., 1941, ApJ, 94, 511 
Chandrasekar S., 1943a, Rev. Mod. Phys., 15, 1 

Chandrasekar S., 1943b, ApJ, 97, 255 

Chandrasekar S., 1943c, ApJ, 97, 263 

Chandrasekar S., 1943d, ApJ, 98, 25 

Chandrasekar S., 1943e, ApJ, 98, 47 

Chandrasekhar S., 1944a, ApJ, 99, 47 

Chandrasekhar S., 1994b, ApJ, 99, 25 

Chandrasekhar S., von Neumann J., 1942, ApJ, 95, 489 



A. Del Popolo: On the nature of gravitational forces 19 

Chandrasekhar S., von Neumann J., 1943, ApJ, 97, 1 (CN43) 

Dominguez-Tenreiro, R., Gomez-Flechoso, M. A., 1998, MNRAS 294, 465 

Eddington, A.S., 1916, M.N.R.A.S., 76, 572 

Elson R., Hut R, Inagaki S., 1987, ARA&A, 25, 565 

Gilbert, I., 1970, Ap.J., 159, 239 

Hernquist, L., 1987, Ap.J. Supp. Ser. 64, 715 

Holtsmark P.J., 1919, Phys. Z., 20, 162 

Kandrup H.E., 1980a, Phys. Rep., 63, n. 1, 1 

Kandrup H.E., 1980b, ApJ, 244, 1039 

Kandrup H.E., 1983, Ap.& S.S., 97, 435 

Kashlinsky A., 1986, ApJ, 306, 374 

Kashlinsky A., 1987, ApJ, 312, 497 

Liddle A.R., Lyth D.H.,1993, Phys. Rep., 231, n 1, 2 

Maoz E., 1993, MNRAS 263, 75 

Peebles P.J.E., 1980, "The large scale structure of the Universe", Priceton University Press, 
Princeton 

Prigogine, I., Severne, G., 1966, Physica, 32, 1234 

Sarazin C, 1988, "X-ray emission from Clusters of Galaxies", (Cambridge: Cambridge Univ. 

Press) 

Seguin, P., Dupraz, C, 1996, A & A, 310, 757 
Strauss M.A., Willick J.A., 1995, Phys. Kept., 261, 271 
White S.D.M., Briel U.G., Henry J.P., 1993, MNRAS, 261, L8 
Wybo M., Dejonghe H., 1996, accepted A & A 295, 347 

Zwart S.F.P., Tout C.A., Lee, H.M., 1998, HiA 11, 622 (Highlights of Astronomy Vol. 11, Kluwer 
Academic Publishers, ed. Johannes Andersen) 



